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Abstract 

Given matrices A and B and vectors a, b, c and d, all with non-negative entries, 
we consider the problem of computing min{c T x : x € Z™ , Ax > a, Bx < b, x < d}. 
We give a bicriteria- approximation algorithm that, given e G (0,1], finds a solution 
of cost 0(m(m)/e 2 ) times optimal, meeting the covering constraints (Ax > a) and 
multiplicity constraints (x < d), and satisfying Bx < (1 + s)b + (3, where (3 is the 
vector of row sums (5i = ^ ■ Bij. Here m denotes the number of rows of A. 

This gives an 0(lnm)-approximation algorithm for CIP — minimum-cost cover- 
ing integer programs with multiplicity constraints, i.e., the special case when there 
are no packing constraints Bx < b. The previous best approximation ratio has been 
0(ln(maxj Aij)) since 1982. CIP contains the set cover problem as a special case, 
so 0(lnm)-approximation is the best possible unless P=NP. 

Key words: covering/packing integer programs, set cover, approximation 
algorithms, multiplicity constraints. 



1 Introduction 

We consider integer covering/packing programs of the following form: 

Given V = (A, B, a, b, c, d) with A G R™ xn , B e R r + Xn , a e M™, b e R r +! and 
c,d G compute OPT = min{c T a; : x G Z", Ax > a, Bx < b, x < d}. 
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The constraints Ax > a, Bx < b, and x < d are called, respectively, covering, 
packing, and multiplicity constraints. 

The width, W, is mmfai/Ay : > 0}. Note that it is easy to reduce any 
instance to an equivalent instance with width W at least 1 — simply change 
each Aij to mm{Aij, a«}. This does not change the set of integer solutions. 

The dilation, a, is the maximum number of covering constraints that any 
variable appears in. 

A p- approximate solution is a solution meeting all constraints and having cost 
at most p times the optimum. A p- approximation algorithm is a polynomial- 
time algorithm that produces only p-approximate solutions. The quantity p is 
called the approximation ratio of the algorithm. 

Perhaps the most well-known problem of the form above is set cover: given a 
collection of sets with costs, choose a minimum-cost collection of sets such that 
every element is in a chosen set. In the corresponding formulation A^ G {0, 1}, 
and cij = 1, for i = l,...,m,j = l,...,n. This problem admits a simple 
(1 + In m)-approximation algorithm [9,12,3], and no o(lnm) -approximation is 
possible in polynomial time, unless P=NP [19]. 

Other special cases include natural generalizations of set cover, including set 
multicover where G Z + and multiset multicover where in addition A^ G Z + . 
[24]. In these problems, multiplicity constraints limit the number of times a 
given set or multiset can be chosen. In facility-location problems (where Xj 
represents the number of facilities opened at a site j), multiplicity constraints 
are used to limit the number of facilities opened at a site. The motivation may 
be capacity limits, security goals, or fault-tolerance (to ensure that when a site 
is breached or damaged, only a limited number of opened facilities should be 
affected) [23,14]. 

We give bicriteria approximation algorithms. For any e G (0, 1], our first algo- 
rithm finds a solution x such that Ax > a, Bx < (1 + e)h + f3,x < |~(1 + e)d] , 
where f3 is the vector of sums of rows of B: j3i = J2j The cost c T x is 
0(1 + ln(l + a) /(We 2 )) times the optimum of the standard linear program- 
ming (LP) relaxation. Note that the standard LP relaxation has an arbitrarily 
large integrality gap if multiplicity constraints are to be respected. Our second 
algorithm finds a solution x of cost 0(1 + ln(l + o)/t 2 ) times the optimum, 
satisfying Ax > a,Bx < (1 + e)b + f3,x < d, thus meeting the multiplicity 
constraints. 
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These algorithms are appropriate for the case when B has small row sums 
(for example, a multiset multicover problem with restrictions such as "from 
the 5 sets si, S2, . . . , S5, only 100 copies can be chosen") and for the CIP 
(covering integer programming) problem, formed by instances without packing 
constraints (no u Bx < b"). CIP is well-studied in its own right. For this 
problem, our second algorithm is an 0(ln(l + a))-approximation algorithm. 
This is the first approximation algorithm for CIP whose approximation ratio 
is logarithmic in the input size. Fig. 1 has a table of known approximation 
algorithms for CIP. 1 

We use here results for another special case — CIP without multiplicity con- 
straints. This problem, which we denote CIPoo, has a long line of research, 
but we use only the following results. Randomized rounding easily yields an 
0(1 + \n(m)/W + \J\n{m) /^^-approximation algorithm, where W, called the 
width of the problem instance, is max{oj/ A^ : A^ > 0}. Srinivasan gives an 
0(l + \n(l + a) /W + \J\n(l + ^/^-approximation algorithm, where a, called 
the dilation of the instance, is the maximum number of constraints that any 
variable occurs in [22,21]. Neither of these algorithms return solutions that 
are suitable for CIP, as the solutions can violate the multiplicity constraints 
by a large factor. 

A preliminary version of this paper appeared in [11]. Other work on covering 
problems includes [4,6,18,15,26,22,21]. See [8] for a survey. 

The outline of this paper is as follows. In Section 2 we present our first main 
algorithm that violates the multiplicity constraints by a (1 + e) factor. In 
Section 3 we discuss the integrality gap of the standard LP formulation and 
present our second main algorithm which meets the multiplicity constraints. 
We conclude in Section 4 with some open questions. 



2 Rounding LP relaxations of CIPoo and CIP 

The approximation ratios in this paper are proven with respect to various 
linear programming relaxations of the problems. Our first main result follows 

1 In the table, H(t) is the harmonic series with t terms. It is well-known that 
H(t) = In £ + 0(1). To give some intuition for the Fisher- Wolsey bound consider for 
example the case where each Cj = 1 and the minimum non-zero entry of A is 1. In 
this case the bound is asymptotically equal to Dobson's. 
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who 


restriction 
on CIP 


cost approximation 
ratio 


multiplicity 
guarantee 


Fisher 

& Wolsey [6] 


none 


l + HPi/fo) 

(3i = maxjJ2iAij/cj 
(3 2 = mm{^\ Aj > 
0} 


x < d 


Dobson [4] 


G Z + 


//(max^E^A,) 


x < d 


Rajagopalan 
& Vazirani [18] 


A^ e {0, 1} 


0(ln(l + a)) 


x < d 


Srinivasan 
& Teo [23] 


Cj = 1 


0{l + ln{m)/{We 2 )) 


x < \{l+e)d] 


Kolliopoulos [10] 


A^ G 
{0,0,} 
for some (f>j 


0(ln(l + a)) 


x < \l2d] 


Srinivasan 
[22,21] 


none 


0(l + \n(l + a)/W) 


x < 0(1 + 
ln(l + a) /W)d 


this paper 


none 


0(l+ln{l+a)/{We 2 )) 


x < \{l+e)d] 


this paper 


none 


0(ln(l + a)/e 2 ) 


x < d 



Fig. 1. Approximation algorithms for the CIP problem, 
min{c T x : x <G 7L\,Ax > a, x < d}. The width W is mm{cii/Aij : Aij > 0}. 
Without loss of generality, W > 1. The dilation a is the maximum number of 
constraints any variable appears in. The algorithms presented in this paper gener- 
alize to allow packing constraints (Bx < b); for the general case the approximate 
solution x satisfies Bx < (1 + e)b + (3 where /3j = J2j Bij- 

from careful consideration of the relation between various forms of the problem 
and their standard relaxations. 

We begin by describing a standard approximation algorithm for CIPoo. Given 
an instance V = (A,a,c) of CIPoo, the standard linear programming (LP) 
relaxation is FOPToo = min{c T a; : x G M™ , Ax > a}. We call feasible solutions 
to this LP fractional solutions to V . In contrast we call actual solutions to V 
integer solutions. 

The value FOPToo can be computed in polynomial time (using linear pro- 
gramming) and is a lower bound on the optimum value opt. The algorithm 
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computes an optimal solution x (of cost FOPToo) to the fractional relaxation, 
then rounds x to an integer solution x using the following randomized rounding 
scheme: 

Lemma 1 (folklore) Given a CIP^ instance V = (A, a, c) and fractional 
solution x, let L — 1 + max{4 ln(2m)/W / , ^4 \n{2m)/W}. With positive prob- 
ability, the following rounding scheme produces an integer solution x of cost 
at most 2L times the cost of x: 

1. Let x' = Lx. 

2. Randomly round x' to x: 

let Xj = \x'j] with probability x'j — [x'A , and Xj = [x'A otherwise. 

The proof is standard and we postpone it until the appendix. In what follows 
the floor (ceiling) of a vector t denotes the vector where the ith coordinate is 
the floor (ceiling) of ij. 

Corollary 2 Given a CIPqo instance V = (A, a, c) and fractional solution x, 
let L — 1 + max{41n(2m)/W / , ^4 \n(2m)/W} . One can compute in polynomial 
time an integer solution x < \Lx] of cost at most 2L times the cost of x. 

The corollary follows because the rounding scheme can be derandomized using 
the method of conditional probabilities [5,16,20]. The rounding scheme above 
has been improved by Srinivasan, who shows the following: 

Theorem 3 ([21]) Given a CIPoo instance V = (A,a,c) and fractional so- 
lution x, let a be the maximum number of constraints in which any variable 
appears. For some L = 1 + 0(ln(l + a) /W + ^/ln(l + a) /W), one can compute 
in polynomial time an integer solution x < \Lx] of cost O(L) times the cost 
ofx. 

Since the optimal fractional solution x can be computed in polynomial time, 
Srinivasan immediately obtains an 0(L)-approximation algorithm for CIPoo. 

2.1 Extending to CIP using 1/ K -granularity 

A natural idea would be to extend the rounding schemes above for dP^ to 
handle CIP problems too. Of course, to do this, we need to figure out how to 
handle the multiplicity constraints. The natural LP relaxation of CIP is 

fopt = min{c T a; : x G K" , Ax>a,x< d}. 
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The first idea would be to compute the optimal fractional solution x, then use 
the rounding scheme from Lemma 1 or Theorem 3 to find an integer solution 
x approximating x. But those rounding schemes return x such that x ~ Lx. 
So, x would violate the multiplicity constraints by a factor of L. But L can 
be as large as fi(mm), and we would prefer to not violate the multiplicity 
constraints so much. 

To work around this, given a CIP V = (A, a, c, d), we do compute an optimal 
fractional solution x, but then, instead of computing an integer solution x that 
approximates x, we first compute a fractional solution x' that is what we call 
(1 / K)- granular — meaning that each coordinate of x' is an integer multiple of 
1/K. We do this for a sufficiently large integer K, so that the (l/X) -granular 
solution x' has x' « (1 + e)x (and satisfies all covering constraints). To get 
the final integer solution x, we round x' up deterministically by rounding each 
coordinate up to its nearest integer. Then x = \x'] < |~(1 + e)x]. A little 
thought shows that this last rounding step increases the cost by at most a 
factor of K, so that the cost of x is 0(K) times the cost of x. 

The next lemma captures the exact tradeoff between granularity and approx- 
imation of the cost (and, implicitly, multiplicity constraints). The lemma is a 
straightforward consequence of the previous results. 

Lemma 4 Fix any integer K > 0. Given a CW^ instance (A, a, c) and frac- 
tional solution x, let a be the maximum number of constraints in which any 
variable appears. For some L = 1 + 0(ln(l + a) / KW + ^/ln(l + a)/KW), one 
can compute in polynomial time a (1/ ' K)- granular solution x" < \Lx] of cost 
O(L) times the cost of x. 



PROOF. Here is the algorithm. The input is V = (A, a, c), x, and K. 



1. Construct CIPqo instance V = (A,Ka,c). Let x' = Kx. 

2. Let x' be the integer solution obtained by applying Theorem 3 to V 
and x' . 

3. Return x" = x'/K. 



Step 2 is well defined as x' is a fractional solution to V'. 

By Theorem 3, x' < \LKx\ is an integer solution to V of cost Q (KL) times 
the cost of x, with L = 1 + (9(ln(l + a)/KW + ^ln(l + a)/KW). 
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Thus (using x" = x'/K), x" < \Lx~\ is a (l/X) -granular solution to V of cost 
0(L) times the cost of x. (We also use here \LKx\/K < \Lx] for integer 
K.) □ 

Note: In Step 2 of the algorithm in the proof, Lemma 1 can be used instead 
of Theorem 3, in which case the 1 + a's in the definition of L (in the lemma) 
are replaced by m's. 

In the remainder of the section, by a (p, £)-bicriteria approximate solution 
for a CIP, we mean an integer solution x that satisfies Ax > a and x < 
\£d], with cost at most p times the optimum FOPT. By a (p, I) -bicriteria 
approximation algorithm, we mean a polynomial-time algorithm that returns 
(p, £)-approximate solutions. 

Our first algorithm works as follows. It first computes a (l/if)-granular so- 
lution x! (where K « ln(l + a) /(We 2 )) approximating the optimal fractional 
solution x. Then it gets an integer solution x by deterministically rounding 
each coordinate of x' up to the nearest integer. It returns x. 

Here is a sketch of the analysis. For this choice of K, x' — (1 + 0(e))x, so 
that x nearly satisfies the multiplicity constraints: x < \(1 + 0(e))x\. Since 
x' meets the covering constraints, so does x. Finally, x' has cost 1 + 0(e) 
times the cost of x, and, crucially, since x' is (l/K)-granular, deterministically 
rounding x' up increases the cost by at most a factor of K. So the final integer 
solution x has cost at most K times the cost of x 1 ', i.e., 0(K) times the cost 
of the original fractional solution x. 

The next lemma gives a detailed statement of the result and its proof. 

Lemma 5 Fix any e G (0, 1]. Given a CIP instance (A, a, c, d) and fractional 
solution x, one can compute in polynomial time an (0(l+\n(l+a)/(We 2 )), 1 + 
e)-bicriteria approximate solution x < \(l + e)x] . 

PROOF. Here is the algorithm. The input is V = (A, a, c, d), x, and e. 

1. Take K = [ln(l + a) /We 2 ] . 

2. Obtain a (l/K)-granular solution x' by applying Lemma 4 to the CIPqo 
instance V' = (A, a, c) with fractional solution x. 

3. Return x = \x'~\ . 
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that x' < \Lx] and that x' has cost O(L) times the cost of x. 

It follows that x < \Lx \ and that x' has cost 0(KL) times the cost of x. (The 
latter because x' is (l/ii')-granular, which implies that the cost of x is at most 
K times the cost of x'.) Since (by the choice of K) L = 1 + 0(e), this implies 
the result. □ 



Remark 1 The result of the lemma is best possible in the following sense. For 
any finite p, a (p, 1) -approximate solution w.r.t. fopt is impossible because of 
the arbitrarily large integrality gap (see Section 3 for an example). It is also 
well-known that the integrality gap for FOPToo is Q(\nm) for the special case 
of set cover where arbitrarily large values for the variables are allowed. Hence 
for any I, a (p, /) -approximate solution for a CIP with p = o(lnm) is also 
impossible. 

Now we can state our first main result — an approximation algorithm for any 
general integer covering/packing problem with multiplicity constraints: 

OPT = min{c T a; : x G Z™ , Ax > a, Bx < b, x < d}. 

The algorithm returns a solution that meets the covering constraints, approxi- 
mately meets the multiplicity constraints (and hence approximately meets the 
packing constraints), and has cost 0(K) times the cost fopt of the fractional 
solution. 

Theorem 6 (first main result) Lets G (0, 1], and an integer covering /packing 
program opt = min{c T x : x G Z™, Ac > a,Bx < b, x < d}, with fractional 
solution x, be given. Let = J^jBij. Then one can compute in polynomial 
time an x G 27? such that 



1. c T x < 0(1 + ln(l + a) /(We 2 )) c T x 



2. Ax > Ax 



> a 



3. x < \(l + e)x] 

4. Bx < (l + e)x + /3 



< \(1 +e)d], and 

< (l+e)b + p. 
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PROOF. Here is the algorithm. The input is V = (A, B, a, b, c, d), x, and e. 



1. Let x be the approximate solution obtained by applying Lemma 5 to 
the CIP instance V' = (A, a, c, d), and fractional solution x. 

2. Return x. 

Properties 1-3 of x follow immediately from Lemma 5. To see that property 4 
holds, note that, from x < \(1 + £)x\ it follows that Xj < (1 + e)xj + 1, which 
implies (Bx)i < (B(l + e)x)i + fa. □ 

The optimal fractional solution x to the LP relaxation can be computed in 
polynomial time, so Theorem 6 immediately implies that the desired approxi- 
mate solution x (having properties 1-4 from the theorem and cost 0(1 + ln(l + 
q;))fopt) can be computed in polynomial time. 

Remark 2 Note that for a CIP problem with max^ dj = 0(1), by taking e = 
1/(2 maXj dj), the above theorem implies that one can find in polynomial time 
an integer solution having cost 0(1 + ln(l + a)/W)FOPT and£j < dj + 1. That 
is, the multiplicity constraints can be met within an additive 1. 

3 Meeting the multiplicity constraints 

Given a fractional solution x, it is not in general possible to find an integer 
solution x meeting the covering and multiplicity constraints exactly and having 
cost 0(ln(l+o;)) times the cost of x. To see this, fix 5 > arbitrarily small, and 
consider the following CIP, which is a simple instance of Minimum Knapsack: 

min{a;2 : x G Z^_, (1 — 5)x 1 + x 2 > 1, X\ < 1}. 

The optimal fractional solution has cost 5, whereas the optimal integer solu- 
tion has cost 1. This example demonstrates that the integrality gap can be 
arbitrarily large if multiplicity constraints are to be respected. 2 However, no- 
tice that the two constraints ((1 — 5)xi + x 2 > 1 and x\ < 1) imply a third: 
x 2 > 5. This third constraint, and the observation that x 2 € Z, imply 5x 2 > 5. 

The constraint u 5x 2 > 5" above is a valid inequality for the CIP, meaning 
that it holds for all feasible integer solutions. Adding a valid inequality to 

2 A similar example appears in [2]. In [18] the integrality gap was erroneously 
claimed to be ^(max^ =1 YliLi Aij). 
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the integer program (IP) does not change the space of solutions or the value 
of the optimal solution. But adding the constraint can strengthen the linear 
programming relaxation by ruling out some fractional solutions, and this can 
give a better bound on opt. For example, adding the constraint to the example 
above, and then solving the LP relaxation with the added constraint, gives a 
lower bound of 1 on opt. 

For the general problem, reasoning as above leads to a class of valid inequalities 
called Knapsack Cover (KC) inequalities. These inequalities generalize valid 
inequalities used for CIP problems with Aij G {0,1} in [1,7,25]. They were 
also used by Carr et al. [2]. 

Our next algorithm begins by finding a fractional solution x to the LP relax- 
ation with a number of KC inequalities added. It then rounds x to an integer 
solution x as follows: for j such that Xj > dj/ (1 + e), it "pins" Xj = dj. (This 
increases the cost by at most 1 + e.) To set the remaining x/s, it rounds the 
corresponding Xj's using the randomized rounding algorithm from (Lemma 1) 
or Srinivasan's algorithm (Theorem 3). Since each non-pinned Xj is at most 
dj/(l + e), this rounding can be done so that Xj is at most dj. 

An astute reader may ask whether this process will work if started with a 
fractional solution x to the LP relaxation without KC inequalities. If so, this 
would yield a faster algorithm. After we describe and analyze the algorithm 
sketched above, we discuss this question. 



3.1 The KC inequalities 

Fix a problem instance V = (A, B,a,b,c,d). For each constraint (Ax)i > a« 
and any subset F of the j's (corresponding to x/s that we imagine pinning), 
define af = max{0, cij — J2jeF Aijdj}. Define also A F = mm{Aij, af } for j G F 
and Af- = for j G" F. In words, af is the residual covering requirement of 
the ith constraint if all variables in F were to be set to their upper bounds, 
and Afj is A^, possibly lowered to ensure the width is at least 1. (In the small 
example above, we knew that, for x<i G Z+, the inequality X2 > 5 held if and 
only if the inequality 5x 2 > S did, so we replaced the former constraint with 
the latter.) The KC inequalities for a set F C iV are A F x > a F . The LP-KC 
relaxation of V is to find x G minimizing c T x subject to Ax > a, Bx < b, 
x < d, and subject to the KC inequalities for all F C N. 

We are not aware of an algorithm that solves this relaxation exactly in poly- 
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nomial time. Carr et al. [2] define the following type of solutions, which are 
adequate for our purpose. For A > 1, call a vector x a X-relaxed solution to 
LP-KC if it has cost at most the fractional optimum of LP-KC and satisfies 
(i) Ax > a, (ii) Bx < b, (iii) x < d and (iv) the KC inequalities for the set 
F\ = {j : Xj > dj/X}. The following theorem follows from the results in [2] 
together with the properties of the ellipsoid method (see, e.g., [13]). 

Theorem 7 ([2,13]) Suppose V = (A, B,a,b,c,d) has rational coefficients. 
For any constant A > 1, a \-relaxed solution to the LP-KC relaxation of V 
can be found in polynomial time. 

For the sake of completeness we sketch the idea behind the theorem. When 
the ellipsoid method queries the separation oracle with a point x, the oracle 
returns a separating hyperplane corresponding either to a constraint of the 
standard LP, or to one that is a valid KC inequality for the set of variables in 
x that are high (in this particular x). In the end, look at the set of hyperplanes 
the separation oracle has passed to the ellipsoid method. That set defines a 
polytope which is a relaxation of the LP-KC polytope. 

The input to our next algorithm is an instance V = (A, B, a, b, c, d) of the 
general problem and an e G (0,1]. The algorithm can also be viewed as a 
reduction of the problem of finding a p-approximate solution to a CIP to 
finding a (p, £)-bicriteria approximate solution for appropriate £. 



1. Set d' := [d\. 

2. Let x be a (1 + e)-relaxed solution to the LP-KC relaxation of V — 
(A,B,a,b, c, d 1 ). 

3. Let F = {j : Xj > + e)}. 

4. Define CIP V = (A',a',c,d") by setting A' := A F , a' := a F , and 
defining fractional solution x' and d" as follows: 

5. For j G F let x) = d] = 0. For j £ F let x) = d'j = Xj. 

6. Find integer solution x' to V 1 by applying Theorem 6 with fractional 
solution x 1 and the given e. 

7. Let Xj = dj for j G F and Xj = x'j for j G" F. Return x. 



Theorem 8 (second main result) Given e G (0,1] , and an integer cover- 
ing/packing program opt = min{c T a; : x G Z™, Ax > a,Bx < b, x < d}, let 
Pi = J^jBij. The algorithm above computes in polynomial time an x E Z" 
such that 
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1. c T x < 0(1 + ln(l + a)/{We 2 ))0VT 

2. Ax > a, 

3. x < d, and 

4. Bx < (l + e)b + p. 

PROOF. Note that the cost of a; is a lower bound on opt. Observe also that 
Step 1 does not change the space of integer solutions. 

First we bound the cost of the solution x' (to the restricted problem V'). Since 
x satisfies the KC inequalities for the specific set F, the definitions of F, A', 
b', and d" ensure that x' is a fractional solution of V 1 . By definition of A F , the 
width of V is at least 1. Thus, the cost of x' is 0(ln(l + a)) times the cost of 
x', which is also 0(ln(l + a)) times the cost of x, and thus 0(ln(l + «)opt). 

Next we bound the cost of the final solution x. The cost of x is at most 1 + e 
times the cost of x, plus the cost of x' . Thus, the cost of x is 0(ln(l + «)opt). 

Next we verify that x does not exceed the multiplicity constraints. This is clear 
for the pinned variables: Xj = dj for j G F . For the other variables (j G" F), 
we have % = x) < \{\ + e)d'^ = [(1 + e)xj\ < \{l + e)dj-/(l + e)l < dj. 

Finally, Bx < (1 + e)b + (3 follows from Bx < b and x < \{1 + e)x] . □ 

Corollary 9 The integrality gap of the LP-KC relaxation for CIP is 0(ln(l + 
a)). 

3.2 Remarks on the necessity of the LP-KC relaxation 

Consider for simplicity that d! = d. The algorithm starts with a (1 + e)-relaxed 
solution x to LP-KC, "pins" Xj = dj for j with Xj > dj/(l + e), then uses 
an existing bicriteria approximation algorithm to set the remaining variables. 
A natural question is whether the KC inequalities are necessary. Would it be 
enough to start with a fractional solution x to the standard LP relaxation of 
the CIP? 

If we do this, the analysis of the algorithm (as it stands) fails because x! may no 
longer be a feasible solution to V . (Indeed, the problem V' may be infeasible 
with d" defined as it is, or even with d" = dj/(l +e). To see this, consider the 



12 



simple example at the start of the section.) This breaks the argument that 
bounds the cost of x. 

Perhaps the first fix that comes to mind is to modify the algorithm to take 
A[j = Aij instead of A\j = Afj for j £ F. But this doesn't work because the 
resulting V' can have width less than 1, worsening the approximation ratio. 

Perhaps the second fix that comes to mind is to modify the algorithm to, say, 
set d" = dj for j (jL F, then solve V from scratch to obtain a (new) optimal 
fractional solution x" . In Step 7, the algorithm would pass that new fractional 
solution x" to Theorem 6 (instead of x') to compute x' . Since the cost of x" is 
still a lower bound on opt, it would seem that we can again bound the cost 
of x as desired. 

The problem with this fix is that the new fractional solution x" can have 
x" > dj/(l + e) for j ^ F. Indeed, it can have x" = dj for j ^ F. Thus, the 
rounded solution x' from Theorem 6 could violate the multiplicity constraints. 

The natural work-around is to augment F by adding any such j to F, then 
start over by returning to step 4 with the new F. But, as this process may 
repeat many times, it is not clear how one might relate the cost of all the 
pinned variables to OPT. 



4 Open questions 

Can one find in polynomial time an integer solution for CIP with an additive 
1 violation of the multiplicity constraints and logarithmic cost guarantee with 
respect to the standard LP optimum (without KC inequalities)? We have 
shown this is possible for the case max^ dj = 0(1). Is there a faster (possibly 
greedy?) 0(lnm)-approximation algorithm for CIP? 
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Appendix 

PROOF, (of Lemma 1) We prove that a; is a 2L-approximate solution with 
positive probability. It suffices to prove that the probability that any of the 
following events happens is less than 1: 

(1) c T x > 2Lc T x, or (2) (3i) (M^W/a* < W. 

Note that E[x] — x' — Lx, so that by linearity of expectation 
E[c T x] = LE[c T x] = L x (fopt^), 

and 

(Vi) E[(Ax)iW/ at] = L(Ax)iW/ai > LW. 

By the Markov bound, the probability of (1) is at most 1/2. 

Note that each Xj can be thought of as a sum of independent random variables 
in [0, 1] (where we consider the fixed part, L x jJ> t° be the sum of \_x'^\ vari- 
ables each taking the value 1 with probability 1). Thus (by the choice of W) 
(Ax)iW/a>i = J2j AijXjW/cii is also a sum of independent random variables in 
[0, 1]. By a standard Chernoff bound [17], 

Pr[ (Ax)iW/ai < (1 - e)LW] < exp(-e 2 LW/2). 

Taking e such that (1 — e)L = 1, for the choice of L in the rounding scheme, 
exp(— e 2 LW/2) < l/2m. Thus, the above bound implies 

Pr[(Ax)iW/ai < W] < l/2m. 

Thus, by the naive union bound, the probability that (1) or (2) occurs is less 
than 1/2 + m/2m = 1. 

We have proven that the randomized rounding procedure returns a 2L-approximate 
solution with positive probability. □ 
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